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ABSTRACT 

We use data from the Multi-Unit Spectroscopic Explorer (MUSE), recently commissioned at the Very Large Telescope (VLT), to 
study the kinematics and stellar population content of NGC 4371, an early-type massive barred galaxy in the core of the Virgo cluster. 
We integrate this study with a detailed structural analysis using imaging data from the Hubble and Spitzer space telescopes, which 
allows us to perform a thorough investigation of the physical properties of the galaxy. We show that the rotationally supported inner 
components in NGC 4371, an inner disc and a nuclear ring - which, according to the predominant scenario, are built with stars formed 
from gas brought to the inner region by the bar - are vastly dominated by stars older than 10 Gyr. Our results thus indicate that the 
formation of the bar occurred at a redshift of about z = (^^ror bars are derived from 100 Monte Carlo realisations). NGC 

4371 thus testifies to the robustness of bars. In addition, the mean stellar age of the fraction of the major disc of the galaxy covered 
by our MUSE data is above 7 Gyr, with a small contribution from younger stars. This suggests that the quenching of star formation 
in NGC 4371, likely due to environmental effects, was already effective at a redshift of about z = r results point out that 
bar-driven secular evolution processes may have an extended impact in the evolution of galaxies, and thus on the properties of galaxies 
as observed today, not necessarily restricted to more recent cosmic epochs. 

Key words. Galaxies: bulges - Galaxies: evolution - Galaxies: formation - Galaxies: kinematics and dynamics - Galaxies: stellar 
content - Galaxies: structure 


1. Introduction 


Evolved disc galaxies are characterised by an extended outer 
disc, dynamically supported by the rotation of its stars around 
the galaxy centre. These discs often present non-axisymmetric 
structures such as spiral arms or bars. The inner region of disc 
galaxies can contain a number of other components: e.g. a dy¬ 
namically hotter spheroid, or disc-like structures, such as inner 
discs and nuclear ring^ Current scenarios attribute external and 
violent formation processes to the former - such as mergers - 
whereas inner disc-like structures are believed to have formed in 
internal and less disturbing processes. These processes are due 
to the impact of the non-axisymmetric components - particularly 
bars - on the otherwise axisymmetric potential of the disc. Those 
merger-built spheroids are sometimes named classical bulges, 
while the internal processes driven by non-axisymmetric com¬ 
ponents are sometimes labelled as part of secular evolution (see 


e.g. [Binney & Merrifield||1998t |Elmegreen||1998t [Athanassoula 


S ee|Buta 

also |Gadotti 


I2013 


]![20T2 


for an extensive review on galaxy morphology; see 
for a list of possible stellar structures in disc galax- 


|20 13 1 l^rmendy |2()T3| |Sellwood|20 1 4| ) . Secular evolution is ex¬ 
pected to be the dominant set of processes in galaxy evolution 
after the peak of star formation activity that occurs between red- 
shifts 1 and 2, i.e. between about 8 and 10 Gyr ago. Before and 
during this intense period in cosmic history, interactions between 
galaxies are believed to be the dominant process through which 
galaxies evolve. After it, in generic terms, for galaxies outside 
dense environments such as clusters or groups, begins a less agi¬ 
tated period when internal secular evolution processes become 
important, particularly those driven by bars formed in galaxy 
discs. 


One of the most important consequences of secular evolu¬ 
tion is the inflow of gas through the bar to the central regions 
of galaxies, which results in the building of inner stellar com¬ 
ponents, a process seen in both theoretical and observational 


work (see e.g. [Combes & Gerin| 1985 [[Athanassoula 1992b[ Sell- 

wood & Wilkinson[ 19931 

[Finer et al.[1995t 

Sakamoto et al.[1999t 

Rautiainen & Salo[2000[ [Regan & Teuben 

2003[|2004l Athanas- 

soula 2005[ [Sheth et al. [2005[ [Wozniak 

12007 IKnapen 20071 

Gadotti|2009t|Kim et al.|2012l|Cole et al. 

2014[ Emsellem et al.[ 
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|2015|[Sormani et al.|2015| and references therein). The inflow of 
gas may originate inner discs and/or nuclear rings - comprised 
by recently formed stars - at or near the inner Lindblad reso¬ 
nance (ILR) associated to the bar, or rather at the region domi- 


nated by the bar X 2 orbits. A number of studies have found 
clear rings associated to bars (from the pioneering works of 

nu- 

Ser- 

sic & Pastoriza[[1965| 1967 

to e.g. [van den Bosch & Emse 

lem 

1998 

[Falcon-Barroso et al. 

2006t[Sarzi et al.|2006l[Kuntsc 

iner 

et al.| 

[2010[[Peletier et al.[[2012| see also Buta & Combes[[l 

996 

and Maoz et al.|2001 ), many of those rings showing intense star 


formation activity when cold gas is available. Other rings just 
evolve passively if most of the funnelled gas has been trans¬ 
formed into stars and there is no additional fuel for more star 
formation. Recent studies show that both the current and past 
star formation activity at the centres of barred galaxies are en- 
hanced when compared to unbarred galaxies (see e.g. |Ellison| 
|et al. [MTU [Coelho & Gadotti||2011| ), although not all barred 


galaxies conform to this picture. [In addition, different studies 
have uncovered diverse results, likely because of differences in 
sample selection (|Cacho et al.|2014t[Zhou et al.|2Q15| ). For ex¬ 


ample, [Cheung et al.| ( |2015 


found “no significant differences in 
le bulges or the gradients of barred 


the stellar populations of t 
versus unbarred quiescent disk galaxies”]. 

Because these inner components are built from material 
brought from the disc (which is already dynamically and struc¬ 
turally mature enough to develop a bar), they are expected to 
show exponential, or near exponential, light profiles, as discs. 
In order to differentiate these inner structures from the clas¬ 
sical concept of bulges residing in the centres of disc galax¬ 
ies, a new nomenclature has developed, leading to terms such 
as pseudobulges, disc-like bulges or discy pseudobulges (|Kor- 


mendy & Kennicutt[[2004l [Athanassoula[[2005l [Drory & Fisher 

2007} [Fisher & Drory[[2010} Gadotti 

20121 ICheung et al.|2013| 

see also Kormendy & Barentine 2010 

[Mendez-Abreu et al. 2014 


ent inner structures in a single galaxy). Much work on these in¬ 
ner, secularly-built components and their relation to classical, 
merger-built bulges is still ongoing, particularly on their place in 
a hierarchical galaxy formation scenario (see [Kormendy 2015 


Sanchez-Blazquez|2015[|Bournaud|2015] |Falc6n-Barroso 2015 

and references therein). 

Placing secular evolution in the hierarchical paradigm of 
galaxy formation is thus an important topic of research in ex- 
tragalactic astrophysics today. Since bars are the major drivers 
of internal secular evolution (and are common in the local uni¬ 
verse, see_e^gjEsOdg£et^I]20^ [Menendez-Delmestre et al. 


2007[[Masters et al.|2011| ), a major step forward would be to set 

the epoch in the evolution of the universe in which bars first ig¬ 
nited secular evolution processes. Ideally, one would like to be 
able to state for how long is the bar in a given barred galaxy 
influencing its host’s evolution. The age of the stellar popula¬ 
tion in a bar is not necessarily a measure of the age of the bar. 
A young bar, i.e. one that has been recently formed, could very 
well be composed of old stars. Thus, different approaches have 
to be developed to address this problem. Progress has been slow 
from an observational viewpoint (given the tremendous difiicul- 
ties faced; see e.g.[Gadotti & de Souza|2005 [2006} [Perez et al 


2007| [2009} [Sanchez-Blazquez et al.[[2Qll [Seidel et al.[[2Q15[ ). 

However, modem simulations at least indicate that bars, once 
formed, are difficult to dissolve (in the absence of major merg¬ 
ers), unless the disc is extremely gas-rich ([Athanassoula et aT 


Bournaud et al.|2005}[Berentzen et al.|2007t[Kraljic et al 


2005 


2012 ) 
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A powerful way to investigate the onset of bar-driven sec¬ 
ular evolution - from a statistical viewpoint - is to study how 
the fraction of barred galaxies evolves with redshift. After the 
pioneering work by [Abraham et al.[ ( [1999| ), most of the latest 
studies have used images from the Advanced Camera for Sur¬ 
veys (ACS), onboard the Hubble Space Telescope (HST), taken 
with filters such as the F814W and F850LP, which correspond 
to broadband 1 and SDSS z, respectively, i.e. as red as possible 
but still within the optical spectrum (see[Jogee et al.|2004[[Sheth 


[et al.|200^ [Cameron et al.|2010[ [Melvin et al.|2014[ ). While the 

results are not entirely compatible between the different studies, 
there seems to be a consensus in that the fraction of disc galax¬ 
ies with bars declines rapidly and monotonically from the local 
universe to z ~ Q.8. [Sheth et al.[pQ12| ) suggest that this is due 
to an evolution in the dynamical state of discs. With time, discs 
grow massive enough and kinematically cold enough to develop 
bars. As alerted by [Sheth et al.[ ( [2008) ), one major difficulty in 
studying the fraction of barred galaxies at z > 0.8 with the red 
ACS filters is that they correspond to rest-frame wavelengths in 
the ultra-violet part of the spectmm, where bars become exceed¬ 
ingly diffi cult to detect (see [Gil de Paz et al.[[2007[ ). To avoid 
this issue, [Simmons et al.[ ( [2014| ) have used images taken with 
the infrared channel of the HST WFC3 (Wide Field Camera 3), 
with filters such as the F160W, allowing them to investigate the 
fraction of barred galaxies at z > 1. Consistently with previous 
results, they found that the fraction of barred galaxies drops to 
z ~ 1, but then stays approximately constant at about 10 per cent 
up to z ~ 2. The images used by Simmons et al. are from the Cos¬ 
mic Assembly Near-Infrared Deep Extragalactic Legacy Survey 
(CANDELS; Grogin et al.|201 1 [Koekemoer et al.|2011[ ). These 
works then indicate that, statistically, bar-driven secular evolu¬ 
tion became important at z < 1. But the work of Simmons et al. 
also suggests that for a number of galaxies these processes have 
started earlier. 


To understand the formation and evolution of galaxies, it is 
also necessary to understand when galaxies form stars and when 
galaxies cease star formation activities. While there are difficul¬ 
ties in understanding how galaxies in the field stop forming stars 
(see e.g. [Martig et al.|2009] ), in dense galaxy environments, such 
as galaxy clusters, the removal of cold gas from galaxies and 
their surroundings is expected to be efficient in quenching star 
formation (see [Balogh et al.[2()04| ). The formation of galaxy clus¬ 
ters is a long process that has initial phases at redshifts larger 
than about 3, and, although massive clusters have been observed 
at redshifts about unity or more ( [Brodwin et al.|2010|[2013t[Pro-[ 
pris et al.|2015t[Ma et al.|2015] ), it is expected that cluster forma¬ 
tion is a process that goes on well below z ~ 1 ( [Kravtsov & Bor-[ 


gani[2012][Planelles et al.[20f4] ). In fact, the nearby Vir go cluster 
is known for not being dynamically fully relaxed yet ( [Binggeli| 
\\999) . One thus expects that the effects of the cluster environ¬ 
ment on the evolution of its constituent galaxies may evolve with 
time. 


In this paper, we explore the capabilities of the recently com- 
missioned MU SE (the Multi-Unit Spectroscopic Explorer, [Ba-[ 
con et al.|2Q10[ ) integral field spectrograph - installed on the Eu¬ 
ropean Southern Observatory Very Large Telescope (VET) atop 
Cerro Paranal, Chile. MUSE is unique in its combination of field 
size, sensitivity, spatial resolution, wavelength range and spec¬ 
tral resolution. We use MUSE to study the inner V squared of 
the nearby galaxy NGC 4371. We use imaging data from the 
Spitzer and Hubble space telescopes, and, with MUSE, we are 
able to study in fine detail the inner components of the galaxy, 
and put constraints to their formation histories. 
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In the next section, we briefly describe the MUSE observa¬ 
tions and the corresponding data reduction. In Sect. we pro¬ 
vide a general descripition of NGC 4371 and its environment. 
We perform in Sect. a detailed study of its structural proper¬ 
ties, as derived using Spitzer images at 3.6yum - from the Spitzer 
Survey of Stellar Structures in Galaxies (S'^G, Sheth et al.||2010 ) 
- and HST images, with their exquisite image quality. This struc¬ 
tural analysis will prove fundamental in our understanding of the 
information extracted using MUSE. The extraction and analy¬ 
sis of kinematical information from the MUSE data cube is de¬ 
scribed in Sect. while in Sect. we describe our study on 
the ages and metallicities of stars in the MUSE field. In Sect. [7] 
we discuss our results in the context of bar-driven secular evolu¬ 
tion and environmental effects in galaxy clusters. Einally, Sect.[^ 
summarises our main conclusions. Throughout this study we use 
a Hubble constant of Hq = 67.8 km s"^ kpc“^ and Qm = 0.308 in 
a universe with fiat topology ( jPlanck Collaboration et al.|2015| ). 


2. Observations and data reduction 

NGC 4371 was observed as a science verification (SV) program 
[60.A-9313(A), PI: D. A. Gadotti] for MUSE during the nights 
25th and 29th of June, 2014. MUSE covers an almost square 
~ 1' X1' field of view with a contiguous sampling of 0.2" x 0.2", 
corresponding to a massive dataset of about 90 000 spectra per 
pointing. A spectral coverage from 4750A to 9300A is achieved 
in the normal instrument setup at a spectral resolution of 2.5A. 
The observations were performed under a seeing with full width 
at half maximum (EWHM) varying from 0.8" to 0.9". Sky trans¬ 
parency was with thin cirrus during the first night and photo¬ 
metric during the second night. Since we do not aim at using 
information sensitive to the flux calibration of the data, any un¬ 
certainty resulting from the thin cirrus in the first night is irrele¬ 
vant here. Nevertheless, the presence of thin cirrus may limit the 
accuracy of the background subtraction. 

We targetted the central ~ 1' x 1' region of NGC 4371 with 
MUSE (see Eig.[T]). The observations were distributed in two 1 h 
observing blocks with a total integration time of 1 h on source. 
As the galaxy fills the MUSE field, we split the exposures per ob¬ 
serving block into 3 exposures of 600 s, and monitored the sky 
background by observing a blank sky field for 180 sec following 
a SKY-OBJ-SKY-OBJ-OBJ-SKY pattern. To be able to reduce 
the effects of bad pixels and fiat-fielding uncertainties we ap¬ 
plied a small dither (< 1") of the field centre and rotated the en¬ 
tire MUSE field by 90° between each object observation. Erames 
were taken to correct the exposures for bias and dark current, to 
fiat-field the exposures, and to perform illumination correction. 
Wavelength calibration is achieved through a set of different arc 
lamp frames, and the exposures were flux-calibrated through the 
observation of a spectrophotometric standard star. Einally, the 
exposures were also finely registered astrometrically. All these 
frames were observed as part of the MUSE standard calibration 
plan, and we applied the calibration frames taken closest in time 
to the execution of our observing blocks. 

The MUSE pipeline (version 0.18.2) was used to reduce the 
dataset. A detailed description of the pipeline will be presented 
in Weilbacher et al. (in prep.), but we briefly outline the reduc¬ 
tion process here. The process is split up into two steps. In the 
first one, each science frame is calibrated separately to take out 
instrumental effects. This consists in subtracting the bias and 
dark current levels, fiat-fielding the images (including illumina¬ 
tion correction), extracting the spectra of all slices, and perform¬ 
ing the wavelength calibration. Those calibrated science frames 


are combined into the final data cube in the second step. It in¬ 
cludes the process of flux calibration, sky subtraction adopting 
a model sky spectrum computed from the sky field, astromet¬ 
ric registration, correction for differential atmospheric refraction, 
and resampling of the data cube based on the drizzle algorithm 
(see [Weilbacher et al.|2012| ). 


3. Introducing NGC 4371: general properties and 
environment 

NGC 4371 is an early-type galaxy that is in the core of the Virgo 
cluster. It is thus nearby (J 16.9 Mpa [Blakeslee et al.|2009 ), 
bright (Mb ^ -19.2 , see HyperLedfl Marprso t.p ~ -21.2, 


Eerrarese et al. 2006), massive (M^ 10^^-^ M©, Gallo et al. 


20101), and classified as a barred SO (lenticular) galaxy ( jButa 
et al.||2015| ). Buta et al. also mention that NGC 4371 has a nu 


clear ring, a barlens, an inner ring (i.e. with a radius close to the 
bar semi-major axis), as well as ansae at the bar ends, and an 
outer lens (see Eig. [^. Ansae are stellar density enhancements 
found often in barred galaxies; they look like the bar handles 
(see jMartinez-Valpuesta et al.|2007| for a recent study on ansae). 
Barlenses have only recently been identified as an independent 
morphological feature. They are believed to be the face-on pro¬ 


jection of box/peanuts (jLaurikainen et al. 20141 [Athanassoula 

et al.||2014 

Athanassoula||2015 see also Eig. 8 in Gonzalez & 

Gadotti|20 

51). In this case then, barlenses are just the inner parts 


of bars that are more extended than the outer parts in the vertical 
direction, as well as along the bar minor axis in the disc plane. 
As we will see below, our spectroscopic analyses do not reach 
far enough from the galaxy centre to include the outer l ens, the 
ansae and the inner ring mentioned by [Buta et al.lpOlS] ). There¬ 
fore, these components will not be addressed here. Our MUSE 
field covers, however, part of the main disc, the bar, the nuclear 
ring and other inner components discussed in detail below, which 
hold important clues to the formation and evolution history of the 
galaxy. 

The top panel of Eig. ^ shows the S'^G mosaic at an im¬ 
age stretch that highlights the faintest components of the galaxy. 
The sensitivity of this image is approximately (Icr) 27 AB mag 
arcsec"^, or about 1 M© pc“^. Yet there is no evidence for a vi¬ 
olent recent interaction. On the other hand, it weakly suggests 
the presence of a very faint polar ring that could have originated 
from the encounter with a satellite galaxy. Therefore, it seems 
reasonable to assume that no dramatic merger has occurred re¬ 
cently in the galaxy. The SDSS colour composite in the middle 
panel shows clearly how uniform is the optical colour distribu¬ 
tion across the galaxy, with no signs of recent star formation. 
The bottom left panel in Eig. shows that our MUSE exposure 
covers almost the whole extent of the bar and a sizeable part of 
the major disc. The bottom right panel highlights a ring with a 
semi-major axis of 10.4", semi-minor axis of 4.3", and po¬ 
sition angle similar to that of the major disc (see also Eigs.j^and 
1^. It has a width of about 2". This is the nuclear ring mentioned 
earlier, and we will call it the 10" ring. 

In Eig. we show the location of NGC 4371 in the Virgo 
cluster, as projected on the plane of the sky. NGC 4371 is a mem¬ 
ber of sub-cluster A ( jGavazzi et al.|1999] ). It is near the bottom 
of the cluster’s potential well, and near one of its brightest galax¬ 
ies, M87, and its extended X-ray halo of hot gas. The distance 
to NGC 4371, m easured from surface br ightness fiuctuations, is 
16.9 + 0.6 Mpc ( jBlakeslee et al.||2009| ), while that to M87 is 
16.4 + 0.5 Mpc (see [Bird et al.|2010| who arrive at this value by 


http://leda.univ-lyonl.fr/ 
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Fig. 1. Top: reduced mosaic of NGC 4371 displayed using his¬ 
togram equalization in order to highlight faint features. The inner red 
box represents the MUSE field of our exposures, with 1' on a side. 
This image suggests the presence of a very faint polar ring that could 
have originated from a satellite galaxy, but the evidence for it is weak 
at best. No sign of a violent recent interaction is seen. Middle: SDSS 
colour composite of NGC 4371, 6.6' on its longer side. Bottom left: 
central region of the S"^G image with isophotal contours overlaid in red 
and a diagram overlaid in black representing our MUSE exposures. The 
MUSE field is the inner trapezoid; the outer polygon area is used to ac¬ 
quire point sources for the slow guiding system. The bar is clearly seen. 
Bottom right: image obtained via unsharp masking of the S'^G image, 
showing clearly the ring near the centre of NGC 4371. This ring has a 
semi-major axis of 10.4" and a semi-minor axis of 4.3". It has a 
width of about 2", and is delineated by the black ellipse. North is up. 
East to the left in all panels. 


averaging different direct measurements). Projected on the plane 
of the sky, the distance between the two galaxies is about 1.5 de¬ 
grees, corresponding to about 0.4 Mpc. The virial radius of the 
cluster is approximately 1.8 Mpc ( [Hoffman et al.|1980|), and the 
heliocentric radial velocity of NGC 4371, 933 km s“^ ( Cappellari 
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Fig. 2. Top: sky map covering the Virgo cluster. The positions of NGC 
4371 and M87 are highlighted, as well as an area of 16 squared de¬ 
grees around the centre of the cluster. Bottom: expanded view of the 
16 squared degrees box highlighted on the top panel. At the distance 
of the cluster, NGC 4371 is about 0.4 Mpc away from M87, which sits 
near the bottom of the cluster’s potential well. (Adapted with permission 
from the Deep Sky Hunter Star Atlas; see http:Hwww.deepskywatch.coml 
- courtesy of Michael Vlasov). 


etal.|201l|i, is close to the mean Virgo velocity of 1035 km s ^ 
( Mould et al.||2000| ). Altogether, these figures indicate that it is 
very probable that NGC 4371 undergoes significant interaction 
with the cluster environment. The results from e.g. jPetitpas & 
jWilsonj p004| ) and [Brown et al.| ( |2011| ) indicate that the galaxy 
contains very little atomic and molecular gas. This could be the 
result of processes induced by the dense environment where the 
galaxy is located. We will come back to this point in Sect. [7] 
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4. Structural analysis 


In order to achieve as thorough as possible an understanding of 
the spatial variations in spectral properties in galaxies, such as 
kinematic measurements and stellar ages and metallicities, it is 
fundamental to study their structural properties. To put it differ¬ 
ently, a very powerful way to investigate the formation histo¬ 
ries of galaxies is to combine spectral information, such as kine¬ 
matics, with photometric/morphological information, such as the 
structural properties of different stellar components in a galaxy. 
As an example, one would use photometric information to un¬ 
derstand where the bulge dominates in a given galaxy and in this 
way be able to isolate - using only spectra from that region - 
the bulge mean stellar metallicity or age from the corresponding 
parameters of th e disc or the galaxy as a whole (see e.g. |Coelh^ 
& Gadotti|2Qll| ). This section is thus mainly dedicated to inves¬ 
tigate in detail the structural properties of NGC 4371 before we 
dive into the spectral properties derived with our MUSE SV data. 

We take advantage of the fact that NGC 4371 was already ob¬ 
served with powerful imagers such as those onboard the Spitzer 
Space Telescope and the HST. We use the former to derive struc¬ 
tural properties of the major disc and bar of the galaxy. We use 
a 3.6yum image (from the Infrared Array Camera, IRAC, chan¬ 
nel 1) obtained as part of the S^G. The pixel size of the S^G 
image is 0.75'', and the P SF (point spread function) FWHM is 

1.8" ( |Kim et ak 2014[ ). At these wavelengths, dust extinc¬ 
tion and emission have little effect in our measures, so that they 
provide a very accurate representation of the bulk of the stellar 
population in the galaxy. [Meidt et al.| ( |2012| ) find that the contri¬ 
bution from PAHs (polycyclic aromatic hydrocarbons) and hot 
dust is between about 5 and 10 per cent of the integrated light at 
3.6yum. The HST image (downloaded from the Hubble Fegacy 
Archive) is used to study the inner structural components, where 
the pixel size and spatial resolution of the Advanced Camera for 
Surveys (ACS) provide a precious sharp view. The pixel size is 
0.05" and the PSF FWHM is 0.09". Although in this case we 
have to use an optical passband, we use an image taken with the 
F850FP hlter, which has the peak of transmission at about 850 
nm, in order to minimise dust effects. This HST image was taken 
as part of the ACS Virgo Cluster Survey ( |C6te et al.|2004 ). 


4.1. Major structures 

The top panels of Fig. show the radial surface brightness pro- 
hle derived from the S"^ image with ellipse hts to the isophotes, 
using the ellipse task in IRAF. The bottom panels show the 
result of an image d ecomposition with BUDDA ( |de Souza et al.| 
|2004[ [Gadotti|2008) . The ht includes m odels for a photom etric 
bulge, a bar and a type II disc (see e.g. |Erwin et ^|20^ [ and 
references therein, for the different types of disc prohle^, and 
it also takes into account the presence of a bright central point 
source. The bulge and the bar are htted using Sersic functions, 
whereas the disc is htted with a double exponential. A signature 
of the ring shown in Fig. can be seen in the surface brightness 
prohle, i.e. the characteristic bump at a galactocentric radius of 
about 10". One also sees a clearly exponential section in the sur¬ 
face brightness prohle, between the ring and the bar, dominating 
to about 20", which thus appears to be an inner disc. We skip 
for now the discussion on whether this structure we are calling 


We note that 


Erwin et ^ ( 2008 1 do not hnd the type II disc break 


in NGC 4371 with a prohle derived from ellipse hts with hxed position 
angle and ellipticity. We also do not hnd the break in a prohle derived 
from a major-axis cut. That is the only difference between the latter and 
the results from our ellipse hts. 




Fig. 3. Surface brightness radial prohle of NGC 4371. The top pan¬ 
els show only the prohle derived from ellipse hts to the S'^G image. 
The bottom panels also show the prohles of the different model com¬ 
ponents obtained with BUDDA: bulge in red, disc in blue, bar in green, 
and central point source in violet. The prohle corresponding to the total 
BUDDA model is in magenta. The horizontal dotted line marks the surface 
brightness level below which the background noise becomes important. 
The panels on the right focus on the inner 90" (shaded area on the left- 
hand panels). The outer exponential disc, as well as the bump caused 
by the bar are clearly seen. We also point out the bump at a radius of 
10" produced by the ring shown in Fig.[^ and a clearly exponential re¬ 
gion that dominates the prohle between the ring and the bar-dominated 
region, that appears to be an inner disc, given its exponential prohle. 
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Fig. 4. Radial prohles of position angle, ellipticity and the Fourier 
coefficient of the ellipses htted to the isophotes of NGC 4371 in the 
S"^G image, as indicated. The right-hand panels focus on the inner 60" 
(shaded area on the left-hand panels). In the ellipticity prohle we indi¬ 
cate features produced by the 10" ring, inner disc, bar and major disc. 



an inner disc is a barlens, but we note that barlenses usually also 
show exponential prohles. We will return to this question in Sect. 

where the analysis of the kinematical properties of NGC 4371 
will shed light on it. The model for the photometric bulge essen¬ 
tially accounts for all light above the inward extrapolation of the 
disc exponential prohle, except the bar. 

The prohles shown in Fig. were also derived from ellipse 
hts. They show the features in the isophotes’ geometric proper¬ 
ties that correspond to the major disc, the bar and the 10" ring. 
Since the ring has a width of only 2", the plateau in ellipticity 
seen between 11" to 17" does seem to be associated with 
the inner disc. Assuming the major disc to be inhnitely thin, and 
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Table 1. Structural properties of the stellar components in NGC 4371 
as derived from the image BUDDA fit. 


Photometric Bulge 


effective surface brightness 

Me,h 

19.1 

effective radius 

re,h 

11.6 

Sersic index 

nh 

2.0 

ellipticity 

Oh 

0.26 

position angle 

PAb 

91 


Disc 


inner central surface brightness 

/^0,inn 

20.6 

inner scalelength 

^inn 

49.7 

ellipticity 

Od 

0.48 

position angle 

PAd 

92 

break radius 

^br 

156 

outer central surface brightness 

/^0,out 

19.5 

outer scalelength 

^out 

36.5 


Bar 


semi-major axis 

^bar 

34.8 

Sersic index 

^bar 

0.2 

ellipticity 

^bar 

0.51 

position angle 

PAbar 

159 

boxiness 

C 

2.7 

Luminosity/Mass fractions 

photometric bulge/total 

B/T 

0.427 

inner disc/total 

Dinn/T 

0.437 

outer disc/total 

Dout/T 

0.056 

bar/total 

Bar/T 

0.077 

central point source/total 

PS/T 

0.003 


Notes. Luminosity parameters are in units of 3.6pm AB mag arcsec"^. 
Spatial measurements are in units of arcseconds. Position angles are in 
degrees from North Eastwards. Note that the fit also includes a central 
point source. 


taking into account the global maximum in the ellipticity profile 
- which happens to be at the region where the disc dominates - 
one can derive the galaxy inclination with respect to the plane 
of the sky to be 60°. The main structural parameters obtained 
from the fit to the image can be found in Table 




Fig. 5. Surface brightness radial profiles of NGC 4371 derived from 
ellipse fits to the HST image. The right-hand panels focus on the inner 
30" (shaded area on the left-hand panels). Bottom panels have the radial 
axis displayed on a log scale. The vertical dashed lines mark the posi¬ 
tions of the nuclear features seen in the position angle and ellipticity 
radial profiles derived from the HST image (see Fig.[^. 



r (arcsec) r (arcsec) 


Fig. 6. Radial profiles of position angle, ellipticity and the ba Fourier 
coefficient of the ellipses fitted to the isophotes of NGC 4371 in the 
HST image, as indicated. The right-hand panels focus on the inner 5" 
(shaded area on the left-hand panels). The positions of nuclear features 
seen in the position angle and ellipticity profiles are marked by vertical 
dashed lines. 


4.2. Inner structures at high resolution 

The fit to th e HST image was done with GALFIT ( [Peng et ak] 
|2002|[20T0| . This is because GALFIT easily allows an arbitrary 
number of components to be included in the model, and this 
proved to be necessary when we studied the HST image, in order 
to account for the different structural components in the galaxy 
central region. It also allows a very careful modelling of the PSF, 
by using a PSF image derived directly from the ACS image, i.e. 
without the need to assume that the PSF follows any particu¬ 
lar function. This will be important here, as below we will look 
at structural components at the very centre of NGC 4371. The 
structural parameters found for the major disc and the bar in the 
fit of the S^G image were kept fixed in the fit of the HST image, 
of which we focus on the central components. 

Figures!^ andshow results from ellipse fits to the HST im¬ 
age. Given the improved spatial resolution with respect to the 
S'^G image, we can now focus on the galaxy central region. Fea¬ 
tures in the ellipticity and position angle radial profiles can be 
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seen at 0.18" and 0.4" from the centre, in the right-hand panels 
of Fig. This figure also shows a peak in the ellipticity profile 
at about 3". 

In Fig. [7] we take a close look at the stunning features shown 
with HST imaging near the centre of the galaxy. This figure 
shows a colour composite using F850LP and F475W HST im¬ 
ages. The bottom panel zooms in the upper panel, focussing on 
the very central features. The 10" ring is again clearly seen in 
the upper panel. One can also see a few other features in the 
bottom panel. From inside out, starting from the very centre, 
a bright yellow-reddish, slightly less elongated structure with a 
semi-major axis of 0.15". Further out, a dusty ring with semi¬ 
major axis of 0.6". Both these features are at positions roughly 
coincident with the features seen in the position angle and ellip¬ 
ticity profiles (at 0.18" and 0.4", respectively, see Fig. [^. Just 
at the outer rim of the dusty ring, at 0.8", one can also see 
a bluish ring. These dusty disc-like structures were also noted 
by jFerrarese et al.| ( |2006| ) and jComeron et al.| ( |2010| ). We will 
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Fig. 7. Colour composites produced from HST images (F850LP and 
F475W) at the Hubble Legacy Archive website. Top: a clear view of 
the 10" ring. Bottom: the very inner region of NGC 4371. The blue 
point source on the North is ^ 4" from the centre. The bottom panel 
clearly shows distinct features in the very centre, all elongated, with 
major axes positioned along the East-West direction. From inside out, 
a bright yellow-reddish, slightly less elongated structure with a semi¬ 
major axis of ^ 0.15". Further out, a dusty ring with semi-major axis 
of 0.6", and just at the outer rim of it, at ^ 0.8", a bluish ring. These 
dust y disc-like structures were also noted by |Ferrarese et al.| ( |2006| see 


also 


Comeron et al.|2010] ). North is up. East to the left. 


call these structures (except for the less elongated, central com¬ 
ponent within a radius of 0.15") collectively as the nucleus, 
and they are included as a single component in the HST image 
fit with GALFIT. 

Circumscribing all these components there is a fainter elon¬ 
gated structure. It can also be seen in the unsharp mask shown in 
Fig.[^ It has a semi-major axis of 3.2" and a semi-minor axis 
of 1.4". It has a width of about 2". The elongated dark blob 
seen at the centre in Fig. [^corresponds to the structures seen in 
Fig.j^within 0.8". 

The fit to the HST image thus includes models for the ma¬ 
jor disc, the bar, the inner disc, the nucleus, and again accounts 
for a central point source, i.e. the less elongated, central compo¬ 
nent within a radius of 0.15". Note that no extra component is 
necessary to account for the presence of a bulge. Also, the 10" 
ring and the 3" disc-like structure at the centre are not included 
in the fit. The discs are fitted with single exponential functions, 
whereas the bar is fitted with a Ferrer function, and the nucleus 
with a Sersic function. The fit extends to 118", before the disc 
break seen with the S^G image, and thus a double exponential 



Fig. 8. The central region of the F850LP HST image after unsharp 
masking. The 10" ring is clearly visible, as well as another elongated, 
disc-like structure with semi-major axis of ^ 3.2" and semi-minor axis 
of 1.4", as delineated by the red ellipse. It has a width of about 2". 
The structures seen at the very centre of the galaxy in the bottom panel 
of Fig. [7] (within ^ 0.8") are within the elongated dark blob at the centre 
of this figure. 

Table 2. Structural properties of the inner stellar components in NGC 
4371 as derived from the HST image GALFIT fit. 


Inner Disc 


central surface brightness 

/^0,id 

17.5 

scalelength 

^id 

5.0 

ellipticity 

^id 

0.40 

position angle 

PAid 

90 


Nucleus 


effective surface brightness 

l^e,n 

17.8 

effective radius 

re,n 

1.0 

Sersic index 

rin 

1.1 

ellipticity 

Gi 

0.30 

position angle 

PA„ 

90 


Notes. Luminosity parameters are in units of (F850LP) AB mag 
arcsec“^. Spatial measurements are in units of arcseconds. Position an¬ 
gles are in degrees from North Eastwards. Note that the fit also includes 
a central point source. In addition, it also includes the major disc and 
the bar as found with the S'^G BUDDA fit, i.e. the parameters for these 
components in that fit were kept fixed here. 


function for the major disc is not necessary. Results from this fit 
can be found in Table [2| 

In the top right panel of Fig. [^ one can see that the inner 
~ 4" can be described as a pair of two different exponentials. 
The outer one starts at about 1" and extends to about 4", and thus 
seems to be associated to the 3" disc-like component. Within 
1" there is another exponential section, therefore associated to 
the nucleus component. This is consistent with the Sersic index 
found for the nucleus, i.e. = 1.1 (see Table [^, or essentially 
an exponential function. The central point source also makes its 
mark in the bottom right panel of Fig. [^ One sees that the cur¬ 
vature of the surface brightness radial profile changes right at 
the position of the inner vertical dashed line, which marks the 
nuclear feature seen in the position angle and ellipticity profiles 
(Fig.§. 

We now have a good understanding on the stellar compo¬ 
nents hosted by NGC 4371, in the context of their structural 
properties, such as geometrical properties. One interesting point 
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already is the lack of any signature of a massive component that 
could be a classical, merger-built bulge. The components fitted 
by the photometric bulge model in the decomposition of the 
image are now understood to be the inner disc, the 10" ring, 
the 3" disc-like structure and the nucleus. All these components 
appear to be disc-related and not the outcome of a catastrophic 
event. This is not to say that they cannot be formed from ac¬ 
cretion of external material, but a simpler explanation for their 
existence does not involve major mergers - internal processes, 
such as evolution driven by disc instabilities such as the bar, are 
a more likely explanation. 

Erwin et aL] ( |2014| ) argue that NGC 4371 hosts a small clas¬ 


sical, merger-built bulge. This claim is partly based on the evi 
dence that the mean ellipticity in the inner 5" is smaller than that 
of the 10" ring and major disc, which indicates a rounder com¬ 
ponent. Consistently with this picture, Erwin et al. also show, 
using SINEONI data, that there is a significant increase in the 
importance of random stellar motion in this region. On the other 
hand, Eigs.andshow that the 3" component is rather disc¬ 
like, as are the components making up what we are calling the 
galaxy nucleus, such as the dusty and the bluish ring-like struc¬ 
tures clearly seen at the bottom panel of Eig.[^ In addition, our 
HST image fit yields a Sersic index of 1.1 for the nucleus com¬ 
ponent, indicating a disc-like structure (see Table typical un¬ 
certainties in the measurement of Sersic indices are about 20%), 
built from the major disc, since minor mergers are known to re¬ 
sult Sersic indic es significantly larger th an 1 (see|Aguerri et al 


2001 


but see also [Querej eta et aL]|2015| and |Christensen et al. 


2014). We cannot say much about the point source component 


within a radius of about 0.15", however. Physically, this corre¬ 
sponds to about 12 pc, which means it could be a nuclear stel¬ 
lar cluster (see |Boker et al.|2()04] ). We will return to discussing 
whether NGC 4371 hosts a small classical bulge in Sect.j^below, 
where again our kinematical measurements will come to help. 

Eurthermore, from the results of our fits with BUDDA and 
GALFIT we also created what we call structural maps. These are 
simply masks, one for each component fitted, as well as for the 
3" disc-like component and the 10" ring, that indicate in the 
MUSE field where each component dominates the emission of 
light from the galaxy, as seen with the S^G and HST images. 
Thus, each spatial element is assigned to a single structural com¬ 
ponent, namely the one that dominates the light at that spatial el¬ 
ement, above the other components that might as well contribute 
to part of the light in that spatial element. We therefore con¬ 
structed structural maps for the central point source, the nucleus, 
the 3" disc-like component, the 10" ring, the inner disc, the bar 
and the major disc. These maps will be used below in order to 
enhance our understanding and interpretation of the results from 
the MUSE data cube regarding kinematics and stellar content. 


5. Stellar kinematics 

In this section, we first explain the steps taken to extract the kine¬ 
matics information on NGC 4371 from the MUSE data cube, 
and then show the results thus obtained. We also connect these 
results with those from the previous section. 


5.1. Extraction of the kinematic maps 

A high quality analysis of the stellar kinematics requires a min¬ 
imum signal-to-noise ratio (S/N) in the spectra. To derive the 
radial velocity and the velocity dispersion that best represent the 
line of sight velocity distribution (LOSVD) at a given position 
in the galaxy, the corresponding minimum S/N is not a limiting 
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Fig. 9. Sample of spectra extracted from our MUSE data cube. Left- 
hand panels show the region around the Mgb feature, while right-hand 
panels show the region around around the Call triplet. Top panels show 
the spectrum extracted from a region in the Northern part of the bar, 
25" from the centre, with an aperture of 1.5", or 6 spaxels in diameter. 
Middle panels show the spectrum obtained from the East side of the 10" 
ring with an aperture of 1.0". Finally, bottom panels show the spectrum 
from the central spaxel. 


condition. However, the derivation of the higher order Gauss- 
Hermite moments and h 4 of the LOSVD requires high S/N 
(see [van der Marel & Eranx|1993[ see also Gadotti & de Souza 


2005). Using a Voronoi binning scheme ( [Cappellari & Copm 


2003), we binned our data spatially to achieve a minimum S/N 


of 100 at each spatial element. Central pixels surpass this limit 
by far and therefore remain unbinned. In order to avoid contam¬ 
ination at each spatial element by spaxels with very weak signal, 
we also implemented a minimum S/N threshold of 3 for a 
spaxel to be considered. Therefore, the edges of the maps shown 
below from our MUSE data cube are not straight lines, since of¬ 
ten pixels at these edges do not reach the minimum threshold in 
S/N. 

In order to illustrate the quality of the spectra used here, we 
show in Eig.j^a sample of spectra extracted from our MUSE data 
cube, corresponding to different parts of the galaxy. As men¬ 
tioned previously, our spectra cover the range from 4750A to 
9300A but for clarity we show here only two spectral regions: 
one around the Mgb feature and other around the Call triplet. 
While the central spectrum was extracted from a single spaxel, 
the spectrum from the East side of the 10" ring is a median com¬ 
bination of spectra within a 1" aperture, i.e. 4 spaxels in diame¬ 
ter. Likewise, the spectrum from the Northern part of the bar, at 
25" from the centre, corresponds to an aperture of 1.5". These 
apertures correspond to the typical Voronoi spatial bins we used 
at these regions. One sees that this spatial binning is efficient in 
keeping the S/N at a suitable level throughout the field. 

We used the pPXE - penalized pixel fitting - code developed 
by [Cappellari & Emsellem| ( [2004| ) to extract the stellar kinemat¬ 
ics, including higher order moments. The result is a LOSVD de¬ 
scribed by a Gauss-Hermite parametrisation ( jGerhardj 1993 [ [van 


der Marel & Eranx|1993| ), with a measure of the radial velocity 


(v), velocity dispersion (cr) and higher order Gauss-Hermite mo¬ 
ments (hs and h 4 ). The hs and h 4 moments indicate, respectively, 
asymmetric and symmetric deviations from a LOSVD that is a 
pure Gaussian. In other words, hs is basically a measure of the 
skewness of the LOSVD, whereas h 4 is essentially a parame¬ 
ter that quantifies the kurtosis of the LOSVD. They can be used 
to examine the orbital structure of the stellar system in ques- 







































Gadotti et aL: The dawning of secular evolution 


tion. For instance, circular motion results in values that are 
anti-correlated with v. And high values of h 4 suggest the super¬ 
position of structures with different LOSVDs (see e.g. [Bender | 
jet al.|1994| and references therein). These values are determined 
through a fitting process where the routine combines a number 
of template spectra from a previously defined stellar spectra li¬ 
brary to fit the galaxy spectrum of each element. Here we used a 
subset of the MILES single stellar population (SSP) model spec¬ 
tra ( jVazdekis et al.||2010|) with a mean r esolution with FWHM 


of ^ 2.5A ( Falcon-Barroso et al. 201 1| to match the spectral 


resolution in MUSE. The resolution of the model spectra is ad¬ 
justed to the spectral resolution of the data before the fitting 
process. We cover the following range of ages and metallicities 
with the SSP template spectra: 0.1 Gyr to 17.8 Gyr, and -0.40 < 
[Z/H] < -1-0.22, respectively. Throughout this work we assume a 
Kroupa initial mass function (IMF, |Kroupa|200 f ) . A wide wave¬ 
length range was employed in the fits, namely from 4750A to 
8800A, which improves the S/N, as compared to narrow wave¬ 
length ranges. 

Our MUSE spectra virtually reveal no emission lines, except 
for a weak [Nil] emission at 6583A in the inner 1", where the nu¬ 
cleus and the central point source dominate. The [Nil] equivalent 
width in the integrated spectrum obtained by combing all spectra 
within 1" is 0.5A, but Ha is seen only in absorption. We will 
come back to this point in Sect. but we point out here that the 
lack of emission lines means that the results from our pPXF fits 
are straightforward and not complicated in this respect. 


5.2. The kinematics of the different structural components in 
NGC 4371 

5.2.1. Kinematic maps 

Figure [T^ shows stellar kinematics maps of radial velocity v, ve¬ 
locity dispersion cr, and Gauss-Hermite moments hs and h 4 . In¬ 
tensity contours equally spaced in steps of about 0.5 magnitudes 
are indicated. NGC 4371 was also observed with SAURON 
(see [Bacon et al.||2001[|de Zeeuw et al.||2002| ), as part of the 
ATLAS 3D survey ^Cappellari et al.[[20li[ ), covering the inner 
30x40 arcsec. Our MUSE data reveal a multitude of details not 
visible in that former data set. It is easy to contemplate the major 
improvement brought about with MUSE by the same PI and the 
MUSE consortium a decade later ( [Bacon et al.|2010[ ). 

The stellar velocity field exhibits an overall regular rotation 
pattern, with enhanced central regions hinting towards distinct 
kinematic components in those regions. This is confirmed in the 
other maps. The line of nodes is almost aligned perpendicular 
to the bar major axis and only shows a very slight twist in the 
region covered by the MUSE field, namely around 5" off the 
centre and in the bottom part of the field. The overall velocity 
seems to further rise towards the edges of the field, suggesting 
that we do not reach the peak in the rotation curve. 

The stellar velocity dispersion is generally elevated within an 
almost circular region of radius of about 20", with several well- 
defined regions presenting varied patterns in terms of velocity 
dispersion. There are two large regions of higher velocity disper¬ 
sion above and below the centre along the bar major axis, and 
other two large regions of lower velocity dispersion along the 
East-West direction. The very centre does not exhibit the high¬ 
est velocity dispersion, but there seems to be four concentrated 
regions around the centre (and about a few arcseconds from it) 
with higher velocity dispersion, two of such regions aligned with 
the East-West direction, the other two aligned with the North- 
South direction. 


The hs parameter strongly anti-correlates with the radial ve¬ 
locity, in particular in the regions of high absolute radial velocity 
and low velocity dispersion, both in the inner regions, as well as 
at the edges of the field. The two regions showing the highest 
absolute values of ha, at about 10" from the centre both at East 
and West, also show the highest h 4 values. There are also two 
well-defined, concentrated regions showing high absolute val¬ 
ues of ha at a radius of about 1-3 arcsec again both at East and 
West. These regions show counterparts with high radial veloci¬ 
ties in the radial velocity map. The h 4 map also shows a region 
of high values that can be described as an elongated ring with 
semi-major axis reaching out to about 20" along the East-West 
direction. A similar elongated region around the centre, extend¬ 
ing to ~ 5" also shows high h 4 values. Two of the four circumnu- 
clear regions with high velocity dispersion show a hint of having 
relatively low h 4 values. These are the ones along the East-West 
direction. A similar relative drop in h 4 for the other two regions 
along the North-South direction is not so discernible. The large- 
scale structure in the h 4 map also displays asymmetric features 
on its edges. 

Another important kinematical parameter to be considered 
here is Ar. This parameter was introduced by [Emsellem et al. 
( [2007 1 ) as a proxy to quantify the observed projected stellar an¬ 
gular momentum per unit mass. It is defined as: 


Ar = 


{R\v\) 


(r Vv^ -I- 


( 1 ) 


where R is the galactocentric radius. Eigure pH] shows the cor¬ 
responding map. This map shows more clearly now the features 
that produce the pinching along the East-West direction in the 
inner 2-3 arcsec. These features were already noticed above in 
the maps of radial velocity and hs, and they correspond to rel¬ 
atively high angular momentum. This map also shows the high 
angular momentum component at 10-20 arcsec from the centre 
along the East-West direction. 

5.2.2. The bar, box/peanut and barlens 

Now, how do all the features seen in these maps relate to the 
structural components discussed in Sect. ^ The high rotational 
support in the region dominated by the major disc is no surprise. 
The bar shows a remarkable region with high velocity disper¬ 
sion. Because the galaxy is quite inclined (i ^ 60°), these mea¬ 
sures of cr have components of velocity dispersion both at the 
plane of the disc and at the direction perpendicular to it. Could 
these enhancements in cr be caused by a box/peanut? To try 
answering this question we compare our maps with those re¬ 
cently published by [Iannuzzi & Athanassoula] ( [2015[ ). lannuzzi 
& Athanassoula have produced maps similar to ours from simu¬ 
lations of barred galaxies at different evolutionary stages, before 
and after the formation of the box/peanut in the bar. Their figures 
29 and 31 show their simulation GTRlOl at a very similar pro¬ 
jection as how we see NGC 4371. Their velocity dispersion maps 
show enhancements along the bar major axis that very much re¬ 
semble our results for NGC 4371, particularly so after the forma¬ 
tion of the box/peanut. Another signature of box/peanuts pointed 
out bypebattista et al.[ ( [2005] ) , and corroborated by lannuzzi &[ 
[Athanassoula[ ( |2015[ see also [Mendez-Abreu et al. |^08[ ), are 
minima in h 4 seen along the bar major axis at each side of the 
centre. Our h 4 measures along the bar major axis suggest the 
presence of such minima, particularly at the bottom of the field. 
Nevertheless, we admit that this evidence is only suggestive. [Er-[ 
win & Debattista[ ( [2013[ ) do not find evidence of a box/peanut 
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Fig. 10. Radial velocity, velocity dispersion, and h 4 maps for the stellar component in NGC 4371, as indicated. The colour bars on the side 
of each panel indicate the plotted range of the parameter measured. For radial velocity and velocity dispersion these are given in km s"^ The 
isophotes shown are derived from the MUSE cube reconstructed intensities and are equally spaced in steps of about 0.5 magnitudes. North is up. 
East to the left. The nearly vertical dashed line indicates the position of the bar major axis. The horizontal elongated markings indicate the outer 
boundary of the region dominated by the inner disc, at ^ 20" - as derived through the ellipse fits - and the position of the 10" ring. 
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Fig. 11. Map of Ar for NGC 4371. The isophotes shown are derived from the MUSE cube reconstructed intensities and are equally spaced in steps 
of about 0.5 magnitudes. The right-hand panel focusses on the inner 10". North is up. East to the left. The nearly vertical dashed line indicates the 
position of the bar major axis. The horizontal elongated markings indicate the outer boundary of the region dominated by the inner disc, at ^ 20" 
- as derived through the ellipse fits - and the position of the 10" ring. 
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in NGC 4371 but they cannot rule out its existence. Overall, it 
seems plausible that NGC 4371 has a box/peanut. Such struc¬ 
tures are known to reduce the strength of bars in inducing secular 
evolution processes and might significantly infiuence the evolu¬ 
tionary processes in a galaxy, such as gas fiows (see |Fragkoudi| 
let al.|2015D . 


As mentioned above, recent work have put forward argu¬ 
ments suggesting that box/peanuts have a counterpart in the 


plane of the disc, which are barlenses ( 

Laurikainen et al.|2014, 

Athanassoula et al.||2014l jAthanassou 

la||2015l jLaurikainen &| 


Salo|2015| ). We do not find evidence for a barlens in NGC 4371 


in the kinematic maps. In the structural analysis done in Sect. 

the only structural component found that could be a barlens 
is what we call the inner disc, which dominates the emission of 
light from the galaxy from a radius of about 10" to about 20". 
Figures [T^ and [TT] show clearly that this component is strongly 
rotationally supported. It is not clear yet whether this would be 
consistent with a barlens, but it is certainly consistent with an 
inner disc indeed. Therefore, if NGC 4371 has a box/peanut, it 
could be that it is not strong enough to produce an evident sig¬ 
nature of its disc plane counterpart, the barlens, in our data. In 
addition, projection effects and/or the presence of the other con¬ 
spicuous inner components in the same region - such as the inner 
disc, in fact - could mask out signatures of the barlens. 


5.2.3. The inner disc, 10" ring and inner structures 


The 10" ring is just at the inner edge of the inner disc and cannot 
be distinguished from it in Figs.[T^and[TT] except perhaps in the 
h 4 map. The fact that this separation is not straightforward hints 
at a connected formation history. At a radius of about 10" one 
sees a peak in h 4 , higher than the values seen between about 
10" and 20", i.e. through the whole extent of the inner disc. The 
peak in h 4 can be explained if the inner disc and the 10" ring 
have different LOSVDs. 

The 3" disc-like component can also be seen in the kine¬ 
matic maps. It is a fast rotating central component, seen as the 
two concentrated peaks, near the centre, of absolute radial veloc¬ 
ity and hs, and Ar. The major disc, the inner disc, the 10" ring 
and the 3" disc-like structure all show a strong anti-correlation 
between radial velocity and h 3 , consistent with circular motion. 
Therefore, all such structural components, found through the im¬ 
age decompositions, have their existence corroborated and better 
understood through the kinematic maps. The central point source 
and the components that we have collectively called the galaxy 
nucleus dominate the central 1". 

Now we can come back to the question on whether NGC 
4371 hosts a small classical bulge, within 5" from the centre, as 
suggested by jErwin et al.| ( |20T4| ). Our maps show clearly that 
the 3" disc-like component is strongly supported by rotation and 
circular motion, which again argues against the presence of a 
classical bulge that extends to a few arcseconds. [Unless, some¬ 
how, a component with relatively high angular momentum re¬ 
sults from the bulge formation processes, despite their violent 
nature. Another possibility is that the small classical bulge did 
not originally have significant angular momentum, but that it 
was spun-up by the bar (see jSaha et JT]|2012| ).] Furthermore, 
we note that, to the spatial resolution limit that we can achieve 
with these observations (i.e. < 1"), the central values of velocity 
dispersion are below those along the inner part of the bar major 
axis, i.e. there is no discernible kinematically hotter component 
in our central kinematic bins. A caveat here, however, is that, 
if the light from the central bins is biased towards a relatively 
younger population that is kinematically colder - even though. 


as we will see shortly, the stellar population content in the centre 
is predominantly very old - this could mask an older, kinemati¬ 
cally hotter component. The measurements presented in Fig. 6 of 
jErwin et^ ( 2014| ) were performed in the near-infrared and with 
adaptive optics, using SINFONI at the VLT, and therefore reach 
a spatial resolution much better than our own measurements (i.e. 

0.1"). From a kinematical perspective, we cannot rule out that 
the nucleus was built through mergers, but we remind the reader 
that our fits to the HST image of NGC 4371 result in a nucleus 
with an exponential, disc-like light profile. In addition, there is a 
peak in ellipticity at a radius of 0.4", where the nucleus dom¬ 
inate. There is also very little we can say about the central point 
source here. 


6. Stellar population content 

As in the previous section, here we will also first focus on the 
methodology employed, in this case to study the properties of 
the different stellar populations in NGC 4371, and then show 
the results thus obtained. We will connect these results with the 
results from our structural and kinematical analyses above, and, 
in particular, make use of the structural maps described in the 
last paragraph of Sect. 

6.1. Extraction of the stellar age and stellar metallicity maps 

In this paper, we have chosen to perform a full spectral fitting 
using the code STECKMAP (STellar Content and Kinematics via 
Maximum A Posteriori likelihood, see |Ocvirk et al.||2006a|b| ). 
The reconstruction of the stellar age distribution and the age- 
metallicity relation is non-parametric, i.e. no specific shape for 
the star formation history is assumed. The problem of obtain¬ 
ing this physical information from the observations is an ill- 
conditioned problem ( |Ocvirk et al.||2006a| ), where small fluc¬ 
tuations in the data can result in large variations in the solu¬ 
tion. To deal with this problem, STECKMAP makes a statistical 
regularisation with the prior in such a way that solutions that 
change smoothly are more likely to result. This is the only a 
priori condition. The smoothness parameter can be set by gener¬ 
alised cross-validation of the possible solutions, accounting for 
the level of noise in the data, in order to avoid over-interpretation 
of the results (see |Qcvirk et al.|2006a| for details). 

The function to minimise is defined as: 


= x"{S {X, Z, g)) + Z, g), (2) 

which is a penalised where S represents the synthetic spec¬ 
trum, X the flux distribution, Z the metallicity distribution, and 
g the broadening function. The penalization can be written 
as P^{x,Z,g) = PxP{^) + PzP(^) + PvP(g), where the func¬ 
tion P gives high values for solutions with strong oscillations 
(flux or metallicity changing rapidly with time) and small val¬ 
ues for smoothly varying solutions. Adding the penalization P to 
the objective function is exactly the same as injecting an a pri¬ 
ori probability density distribution to the solution, as /prior(-^) = 
exp(-/ijcP(-^)), where P is a quadratic function of the u nknown. 

We use again the MILES SSP spectral models (|Vazdekis 
et al.||2010| i, spanning an age range from 6.3x10^ to 1.7x10^^ 
years, divided in 30 logarithmic age bins. The metallicity ranges 
from [Z/77] = -1.3 to [Z/77] = -rO.2. 

Although STECKMAP can be used to obtain the broadening of 
the lines, here we use the values obtained in the previous sec¬ 
tion, calculated with pPXE, and thus fix the broadening function 
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Fig. 12. Mass (top panel) and flux (middle panel) fractions as a function 
of age, and the age-metallicity relation (bottom panel), obtained for the 
central spectrum of NGC 4371. Metallicities are plotted only for flux 
fractions larger than 5%. 


during the fit. The reason for this choice comes from the exist¬ 
ing degeneracy between metallicity and cr ( [Koleva et al.||2008t 
[Sanchez-Blazquez et al.|20lf] ), which biases the mean-weighted 
metallicity if both parameters are fitted at the same time. We also 
do not fit the continuum, to avoid dealing with possible fiux cal¬ 
ibration errors or dust extinction. Instead, we multiply the model 
by a smooth non-parametric transmission curve. This curve is 
obtained by interpolating a spline function between 30 nodes 
spread uniformly along the wavelength range. The choice of fit¬ 
ting or not the continuum may have a strong influence in the de¬ 
rived parameters when the S/N of the spectrum is not very high. 
For spectra with S/N >30 per A the results obtained are inde¬ 
pendent on rectifying or not the continuum (|Sanchez-Blazquez| 
let al.|2011D . 


We ran STECKMAP on each of the binned spectra, and ob¬ 
tained star formation histories and age-metallicity relations for 
all of them. Figure [T^ shows an example of a typical STECKMAP 
output, showing the fiux, mass and metallicity distributions as a 
function of age. 


We can obtain mean values of age and metallicity per spec¬ 
trum weighting with either the fiux or the mass. Flux-weighted 
values are biased towards younger populations, as young stars 
are brighter in the considered wavelength range. Mass-weighted 
values may be, on the other hand, much more uncertain, es¬ 
pecially when a large fraction of old, low-luminosity stars are 
present. To better understand how uncertain are our mean stel¬ 
lar age estimates at each spatial resolution element - which will 
be discussed at length below - we have performed 100 Monte 
Carlo realisations using all binned spectra. We find that the er¬ 
ror in luminosity-weighted age is anti-correlated with the mean 
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age in a logarithmic scale. Considering luminosity-weighted es¬ 
timates, for a mean stellar age of 10 Gyr the median of the er¬ 
ror distribution is 0.035 dex, while for a mean stellar age of 
7 Gyr the median error is 0.05 dex. Therefore, for these mean 
ages, the uncertainty in the mean stellar age is ~ 0.8 Gyr. This 
anti-correlation is not present, however, when we consider mass- 
weighted estimates, which have a median error of 0.014 dex. 
One should note that the Monte Carlo realisations only address 
uncertainties of a Poissonian nature. Uncertainties arising from 
the spectral models employed are not included. In addition, our 
age measures may have as well systematic effects towards higher 
values, as we have measures reaching ages slightly older than 
the currently accepted age of the universe, in particular when 
one considers mass-weighted estimates. In an attempt to neu¬ 
tralise this effect, our main conclusions below are based on the 
luminosity-weighted estimates, which are systematically lower 
than the mass-weighted estimates. In addition, we consider the 
lowest age values whenever there is a wide enough distribution 
of mean stellar ages. 

6.2. Spatial distribution of the different stellar populations 

Figureshows the luminosity- and mass-weighted maps of age 
and metallicity for NGC 4371, derived using the relations above. 
The mass-weighted mean stellar age map shows that the vast ma¬ 
jority of stars in the galaxy are old, with ages above 10 Gyr. The 
major disc is seen to show stars that are on average somewhat 
younger than the rest of the galaxy, but still old. Similarly less 
old stars on average appear also just inwards of the West side 
of the 10" ring. Bins containing the oldest stars are seen pref¬ 
erentially in parts of the bar and the inner disc. The luminosity- 
weighted map shows more structure. It shows again less old stars 
in the disc and in parts just inwards of the 10" ring (i.e. in the 
inner part of the inner disc), especially its West side. The bar 
shows a somewhat intermediate age population, with mean ages 
between those of the major disc and the remaining central com¬ 
ponents. Still, even in the luminosity-weighted map, almost all 
bins inwards of the region dominated by the major disc show 
mean stellar ages above 10 Gyr. 

From the luminosity-weighted metallicity map, it can be 
seen that the disc shows a very homogeneous metallicity dis¬ 
tribution. In addition, the outer parts of the bar show a slightly 
elevated metallicity, while the central regions show a remark¬ 
able peak in metallicity. Finally, the 10" ring shows as well 
elevated metallicity, in comparison to the outer components. 
The mass-weighted map shows less features, but also indicates 
a very strong peak in the central metallicity and the elevated 
metallicity of the 10" ring. 

In Fig.[T^ we show luminosity-weighted and mass-weighted 
maps that represent the location of three different stellar popu¬ 
lations in NGC 4371, according to their ages: young (age < 1 
Gyr), intermediate (1 < age < 7 Gyr), and old (age > 8 Gyr). 
They roughly correspond respectively to a formation redshift of 
around z < 0.08, 0.08 < z < 0.8 and z > 1 (although we stress 
again here that the correspondence between mean stellar age and 
formation redshift is not trivial, given the uncertainties in the 
age estimates). The colour-coded scale corresponds to the frac¬ 
tion of the corresponding stellar population (young, intermediate 
age or old) at each spatial element. Note that the scale changes 
for each panel, in order to show the corresponding relevant fea¬ 
tures. The scale ranges are indicated below the colour bars on 
the right-hand side of each panel. One can see that most stars 
in the central arcminute squared of NGC 4371 are older than 8 
Gyr. Figure [T^ shows that basically one only finds stars formed 
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Fig. 13. Maps for NGC 4371 of mean stellar age (left panels) and stellar metallicity (right panels) weighted by both luminosity (top panels) 
and mass (bottom panels). The age maps have units of Giga-years, and metallicities are given in the spectroscopic notation (logarithmic scale 
normalised to the solar value). 


at redshifts z < 0.08 at the outskirts of the field (i.e. in the ma¬ 
jor disc), and to a much lesser extent in the bar. Stars formed 
between about 1 and 7 Gyr ago are found preferentially in the 
major disc and in the portion of the inner disc inwards of the 10" 
ring. To a lesser extent again, some such stars can be found in 
the bar. These younger stars, however, represent an insignificant 
fraction of the stellar content in this region of the galaxy probed 
by our MUSE observations. One important point to make here is 
that the mean stellar age in all these bins in which stars younger 
than 7 Gyr are found is actually above 8 Gyr, as demonstrated 
in Fig. This will be further illustrated below. The outer part 
of the inner disc and the inner 5" show the highest densities of 
stars formed more than 8 Gyr ago. 

There are other powerful ways of visualising the informa¬ 
tion content in such maps that may disclose important clues in 
a more straightforward way. One of such ways is presented in 
Fig. [T^ which shows the normalised distributions of the mean 
stellar age and metallicity of the spatial bins in each structural 
component. These were derived through a combined use of the 
structural maps constructed in Sect. 4|and the age and metallicity 
maps produced in this section (Fig. T^. These distributions cor¬ 
roborate our previous assessments and also reveal further aspects 
of the stellar content in NGC 4371. The luminosity-weighted 


metallicity distributions appear to be bimodal for the major disc 
and bar components only. Their mass-weighted metallicity dis¬ 
tributions show a peak that roughly coincides with the high 
metallicity peak of their luminosity-weighted distributions, but 
they also show a tail to lower metallicities. Another interesting 
feature clearly seen in the luminosity-weighted age distributions 
is that all components closer to the centre than the 10" ring have 
on average an older population of stars than the remaining of 
the galaxy, particularly the 3" disc-like component. The mass- 
weighted distributions, however, show much subtle differences: 
they show that the whole stellar content of the galaxy is almost 
entirely very old. 


Finally, we present in Fig. plots of the mass-weighted 
and luminosity-weighted mean values of stellar age and 
metallicity for all structural components in NGC 4371. Since 
the components are arranged along the x-axis in a way that 
from left to right each component successively dominates at 
larger galactocentric distances, these plots can be thought of 
as stellar age and metallicity radial profiles. Again, this figure 
repeats some of the information already discussed above, but 
this visualisation format allows a more direct comparison of 
the mean values of each component. In addition, it shows in a 


Article number, page 13 of 


19 
























A&A proofs: manuscript no. N4371_v3 


<lGyr lGyr<age<7Gyr >8 Gyr 



-30 -20 -10 0 10 20 30 -30 -20 -10 0 10 20 30 -30 -20 -10 0 10 20 30 

Aa (arcsec) Aa (arcsec) Aa (arcsec) 


Fig. 14. Top panels: luminosity-weighted maps of stellar populations separated by mean stellar age, as indicated (Young; Intermediate; Old). The 
colour-coded scale, indicated by the colour bar at the right-hand side of each panel, corresponds to the fraction of the stellar population (young, 
intermediate age or old) at each spatial element. Bottom panels: corresponding mass-weighted maps. Note that the scale changes for each panel, 
in order to show the corresponding relevant features. The indicated limits of the colour bars, however, show that most of the mass and light are in 
(or come from) stars older than 8 Gyr. 


straightforward fashion how these values vary from the centre 
to the outskirts of the MUSE field. 


The mean stellar age profiles show how homogeneously old 
is the stellar population in the inner arcminute of NGC 4371 
(particularly so in the mass-weighted results), with the indica¬ 
tion that younger stars are preferentially found in the outer part 
of the field, i n the major disc. [Similar results are found by |Seidel| 
et al. ( |2015| ) for NGC 5701, another early-type barred galaxy.] 
A mild drop in the mean stellar age is also seen at the 10" 
ring. The mean stellar metallicity profiles clearly show the pro¬ 
nounced peak in metallicity at the central components. Again, 
the 10" ring produces a bump in the curves, this time indicating 
a higher metal content than what expected from the global trends 
of the profiles, i.e. as compared to a metal content that grows 
smoothly towards the centre. These profiles show essentially that 
beyond the inner few central arc seconds the mean metallicity 
does not vary much, whereas it starts increasing in the central 
arcsecond until it reaches a maximum at the very centre. The 
luminosity-weighted profiles in Fig. show a somewhat dif¬ 
ferent behaviour when compared to the mass-weighted profiles, 
possibly indicating different properties between a more luminous 
and relatively younger stellar population and the old bulk of the 
stellar content. These differences, along with all the results pre¬ 
sented in this section and in Sects. and will be interpreted 
and discussed in the next section. 
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7. Discussion 

7.1. NGC 4371: a fossil record of the oldest bars 


Being it a disc galaxy, and one that hosts a large-scale bar, one 
expects that internal secular evolution processes have played an 
important role in the evolution of NGC 4371. The BUDDA fit (see 
Table [T]) indicates that the bar-to-total luminosity ratio in the 
galaxy (Bar/T = 0.077) is typica l of the values found for other 
massive galaxies. |Gadotti| ( |2011| ) finds that the typical Bar/T is 
about 0.1 in a sa mple of nearly 300 m assive barred galaxies at 
z ~ 0.05 (see also [Weinzirl et al.|2009 ). Interestingly, the lumi¬ 
nosity profile of the bar in NGC 4371 is very fiat. This is in¬ 
dicated by its low Sersic index, namely TZbar = 0.2. |Kim et ak 


( 2015| ) have conjectured that recently formed bars have expo 
nential luminosity profiles, i.e. with a Sersic index 1, but, as 
bars evolve, their luminosity profiles become flatter, which trans- 
lates to low er values of TZbar- The flattest bars in the work of |Kim| 
et al. ( 201^ , with 144 nearby barred galaxies from the S'^G, have 


a Sersic index similar to that of the bar in NGC 4371. This sug¬ 
gests thus that the bar in NGC 4371 formed at the epoch were 
the first long-standing bars were being formed. Being an evolved 
bar, it makes it more likely that a box/peanut has formed, as such 
structures are seen in simulations to form 1 - 2 Gyr after the 
formation of the bar (see e.g. [Martinez-Valpuesta et al.|2006| ). 


As discussed in the Introduction, one of the main processes 
induced by large-scale stellar bars in disc galaxies is the inflow of 
gas towards the central regions. This inflow of gas is halted near 
the bar ILR and originates nuclear rings and inner discs (also 
known as disc-like bulges) with recently formed stars. We have 
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Fig. 15. Normalised distributions of the mean stellar age and metallicity 
per spatial bin, separated by structural component, as indicated. The 
mass-weighted distributions are plotted in red, whereas the luminosity- 
weighted distributions are plotted in blue. 


suggested above that because the inner disc and the 10" ring in 
NGC 4371 are hardly distinguishable in some of the kinematic 



pointsource nucleus 3"disc-like 10"ring innerdisk bar disk 


Fig. 16. Mass-weighted and luminosity-weighted mean values of stellar 
age and metallicity for all structural components in NGC 4371. The 
structural components are arranged along the x-axis such that this axis 
can be thought of as the galactocentric distance, l.e., the point source 
component is at the centre of the galaxy, and to the right-hand side of it 
we arranged the other components that dominate the emission from the 
galaxy at successively larger distances. The error bars are an outlier- 
resistant measure of the dispersion about the mean, analogous to the 
standard deviation. 


maps, they could have a connected formation history. Therefore, 
the most natural explanation for their origin is that of bar-driven 
gas inflows. Figure [T^ shows that there is no spatial bin in the 
inner disc or the 10" ring with a mean stellar age below 10 Gyr. 
Thus, the last significant gas infiow induced by the bar must have 
happened 10 Gyr ago, and then it is not only the case that the bar 
in NGC 4371 is populated by old stars but the bar itself, as a 
structure, is old. In fact, the mean stellar ages of the inner disc 
and 10" ring imply that the bar was already in place at least 10 
Gyr ago, i.e. at a redshift of about 1.8. Considering the results 
from the analysis of the uncertainties in our stellar age estimates, 
described at the end of Sect. |6.1| the formation redshift of the bar 
is then most likely between z = 1.4 and z = 2.3. 

In the Introduction, we mentioned the study of [Simmons 
et al.| ( [MT^ on the fraction of barred galaxies at high redshifts. 
In that work, the highest redshift galaxy found hosting a bar is 
at z = 1.97, similar to the formation redshift we derive for NGC 
4371. The latter is then a fossil record of those bars, such as 
the ones found by Simmons et al., that formed at early cosmic 
epochs. Interestingly, the bar in NGC 4371 thus testifies to the 
robustness of bars. 

Alternatively, the inner disc and the 10" ring coul d have 
formed via a yet unfamiliar mechanism (see discussion in |Com- 
|er6n et al.|2014[ see also |Eliche-Moral et al.|201l] ), which could 
mean that one cannot constrain the formation redshift of the bar 
using their mean stellar ages. One possibility is that they were 
formed by a bar that has dissolved before the formation of the 
current one. Kraljic et al.| ( |2012] ) show simulations in which bars 
that form at redshifts larger than 1 dissolve, and a long-standing 
bar forms afterwards. However, it is still unclear whether such 
primordial bars would be able to form long-standing inner discs 
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and nuclear rings. It is thus unclear if the 10" ring would survive 
the dissolution of the primordial bar and the formation and evo¬ 
lution of the bar currently observed. On the other hand, Comeron] 
et al.| ( |20f0| ) report a fraction of 19 + 4 per cent of nuclear rings 
found in unbarred galaxies. Nevertheless, the 10" ring in NGC 
4371 is similar to the other 78 nuclear rings in the barred galax¬ 
ies studied by [Comeron et al.| ( |2010| ), in that the ratio between 
the ring radius and the bar length in NGC 4371 is very typical, 
consistent with a picture in which the nuclear ring is formed by 
the bar. Overall, it is thus unlikely that the inner disc and 10" 
ring in NGC 4371 were not formed by the bar. 


7.2. The onset of bar-driven secular evolution in cosmic 
history 

The fact that NGC 4371 is a massive galaxy with an old bar is 
consistent with the results from jSheth et al.| ( |2008||2012| ), which 
indicate that the mor e massive disc galax ies form their bars first. 


For M* > IO^^Mq, Sheth et al. (|2008l find that about half of 


the disc galaxies are barred already at z ~ 0.8. One interesting 
outcome of the present study is that we can provide an estimate 
of the formation redshift for the bar in NGC 4371, and, since 
we know the galaxy stellar mass, we can go a step further and 
use this information as a benchmark to set how the bar forma¬ 
tion epoch evolves as a function of galaxy stellar mass. There¬ 
fore, with a bar formation redshift of z 1.8 and a stellar mass 
M* = IO^^-^Mq, NGC 4371 sets this benchmark and establishes 
that more massive galaxies have formed their bars at even higher 
redshifts. 

It is interesting to consider then, that the disc in NGC 4371 
already had - at these relatively high redshifts - the structural 
and dynamical conditions favourable to the development of a 
long-standing bar. The internal bar-mode disc instability can 
happen spontaneously to any disc galaxy under the right con¬ 
ditions - which are, nevertheless, complex and challenging to 
spot (see e.g. discussion in |Athanassoula|2008}|Sanchez-Janssen| 
|& Gadotti||2Q13| ) - but might also be catalysed by the pertur¬ 
bation of a companion galaxy. A number of studies have been 
dedicated to understand the effects of the environment on bar 
formation (e.g. Aguerri et al.|2009t Li et al.|2009[ Barazza et al 


et al. 


20091 [Mendez-Abreu et al.|2010||20121|Skibba et al.|2012t|Lin 


[2014[ ). While it is interesting to consider the reasonable 


possibility that the bar in NGC 4371 is the result of an interac¬ 
tion, since the galaxy is in a dense environment, this is a difficult 
claim to defend. In any case, there is no observational evidence 
to date that bars catalysed by encounters are in any way different 
from bars formed spontaneously. 


7.3. Environmental effects: the lack of gas and recent star 
formation 

The vast majority of the stars in the central region of NGC 4371 
have ages above 8 Gyr (see Fig.[^, and thus were formed at 
redshifts larger than about 1. In addition, there are no signs of a 
significant recent interaction throughout the galaxy. The galaxy 
thus appears to have gone through no significant accretion of 
cold gas to its central region for the past 8 Gyr. 

As mentioned above, the only emission line we find in the 
spectra is a weak [Nil] line confined within the central 1" and 
with equivalent width of 0.5A. The [Nil] flux could originate 
from stellar photoionization powered by young high-mass OB 
stars, or from other physical processes, such as shock heating or 
photoionization with a power-law spectrum (see e.g. discussion 


in [James et~ST][2005[ ). The central supermassive black hole in 
NGC 4371 - with a mass determined via dynamical modelling 
(P. Erwin et al., in prep.) - is currently not accreting ( [Simoes| 
[Lopes et al.[[2007| ), so we will assume that all [Nil] emission 
comes from recent star formation activity. Even so, the measured 
equivalent width indicates a very low star formation activity cur- 
rently at the centre of N GC 4371 (see Fig. 3 in [Kennicutt|1998[ ). 
[Nagar et S^ ( 2005[ ) and [Capetti et al.[pOQ9[ ) find as well at best 
very low radio emission from the nucleus of NGC 4371, also 
consistent with low star formation or AGN activity and a low 
amount of cold gas. As mentioned at the end of Sect. other 
studies have found very little atomic and molecular gas content 
throughout the galaxy. 

Being it a member of the Virgo cluster, known effects at 
such dense environments could have been responsible for the 
lack of significant star formation and gas content in the cen¬ 
tral region of the galaxy in the past 8 Gyr. As we have seen, 
tidal interactions and mergers with other galaxies do not appear 
to have played a major role. Ram pressure stripping, however, 
could have removed gas from the galaxy, via the interaction be- 
twe en the galaxy cold in ter-stellar gas and the hot intra-cluster 
gas ( [Gunn & Gott| 1 972| ) . If the intra-cluster gas density is high 
enough, e.g. near the centre of a cluster, [Quilis et al.[ p000[ ) 
conclude that ram pressure stripping can remove 100% of the 
atomic gas content in massive galaxies in a time scale as fast as 
10^ yr. The hot intra-cluster gas may also remove any available 
gaseous halo around the galaxy that would o therwise fall into the 
galaxy and feed further star formation (see [Larson et al.[[1980[ 
[Bekki et al.|2002[ ). Such process, dubbed starvation or strangu¬ 
lation, is able to destroy such a gaseous halo even away from the 
cluster core, although in this case on a time scale of 1 Gyr. 
Another effect, dubbed harassment, can also remove an exist¬ 
ing gaseous halo via the weak interactio n with relatively nearby 
galaxies in e.g. a dense cluster core (e.g. [Richstone| 1975} [Moore 
[et al.[pV98] ). The time scales for harassment to be effective in 
changing the stellar population content in a disc galaxy are of 
the order of a few Giga-years. In summary, given the position of 
the galaxy in the cluster, it is not unlikely that NGC 4371 has 
undergone some ram pressure stripping and starvation, which 
would explain the mean old ages of its stars (see e.g. [Vollmer 
[et al.|2001t[Kenney et al.|20T^ . 

Whatever the mechanism that made the inner V squared of 
NGC 4371 devoid of fuel to produce new stars, it was completed 
after the formation of the bar and the subsequent gas inflow to 
form the inner disc and the 10" ring. There is no need to invoke 
ram pressure stripping to cease star formation in the area within 
the bar radius, however. The bar is abl e to push gas to the centre 


fast (at a typical time scale of ~ 10^ yr,[Athanassoula 1992a 


Em- 


sellem et al.|2015[ see also [Huang & kauffmann|2015[[GavazH 


et al.|2015[ ), and the subsequent formation of the inner disc and 


nuclear ring and other nuclear disc-like structures might very 
well be able to consume all available gas. The bar also restricts 
the flow of gas coming from outside its radius in the plane of the 
disc (see [Boumaud & Combes [2Q02| ). However, it is necessary to 
prevent that gas external to the galaxy falls in the disc within the 
bar radius from a direction that is not aligned with the disc plane. 
Starvation/strangulation and harassment come naturally to mind. 

The youngest spatial bins in the disc have a mean stellar 
age of about 7 Gyr (see Fig. [ig, corresponding to a formation 
redshift of z about 0.8. This would mean that the intra-cluster 
gas was already exerting pressure at such early times, when the 
formation of the cluster was still at an early stage. One should 
also consider the likely possibility that strong environmental pro¬ 
cesses could have already occurred if NGC 4371 was part of a 
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group before entering the cluster (see |Haines et al.|2015| ). In any 
case, such environmental effects must have been completed at a 
redshift of about 0.8. 

Since the bar has some spatial bins with mean stellar age 
as low as 8 Gyr (corresponding to a formation redshift of about 
unity, Fig. [Tg, remnant gas promoted some star formation in the 
bar itself after the formation of the inner disc and the 10" ring, 
but before the completed quenching at a redshift of about 0.8. 
In this context, it is also interesting to note that the difference 
in mean stellar age between the peak of the distribution for the 
10" ring and the younger spatial bins present there is about 1 
Gyr. One expects that the gas brought to the centre by the bar 
gets compressed at its ILR and form stars, forming the ring, but 
also originating supernova events that will put some gas back in 
the inter-stellar medium. This may result in a second generation 
of stars, which is probably manifested in the lower age tail of 
the mean stellar age distribution of the 10" ring. This newest 
generation of stars could be the reason why the mean stellar age 
of the 10" ring is slightly less than that of the inner disc (see 
Fig. [^. Furthermore, this second generation of stars may also 
be the reason why the mean metallicity of the 10" ring is slightly 
elevated, as compared to that of the inner disc (again, Fig.[T6|). 

Interaction with the intra-cluster medium thus appears to 
have played a major role in the evolution of NGC 4371, quench¬ 
ing star formation. From a structural perspective, however, the 
galaxy does not present signs of interactions with other galax¬ 
ies or mergers (but see |Wang et ^|2015[ who suggested that 
disc-dominated galaxies without classical bulges can be formed 
via major mergers). In fact, the evolution of the galaxy seems 
markedly governed by internal processes, which is intriguing 
considering its position at the core of the Virgo cluster. One pos¬ 
sibility is that the galaxy is coming near the core of the cluster 
for the first time, but to test this is beyond the scope of this study. 
NGC 4371 is, nevertheless, an example of a galaxy whose evolu¬ 
tion is strongly influenced by internal processes despite residing 
in a dense environment. 


8. Summary and conclusions 

We have used VLT/MUSE integral field spectroscopy data, com¬ 
bined with imaging data from the HST/ACS and Spitzer/IRAC 
imaging cameras, in order to investigate the kinematics, stellar 
population content and structural properties of NGC 4371, an 
early-type massive barred galaxy in the core of the Virgo clus¬ 
ter. The MUSE spectral coverage and resolution, allied to its fine 
spatial resolution, allows us to make robust statements about the 
stellar population content in the inner 1' squared of the galaxy. 
We summarise our main conclusions as follows: 

1. We find that the mean stellar age in the inner disc and 10" 

ring close to the centre of the galaxy is above 10 Gyr, with 
a scattered and negligible contribution from younger stars. 
With the safe assumption that the stars in this ring were 
formed by gas brought to the centre - via bar-induced secular 
evolution processes - and compressed at or near the bar ILR, 
this result indicates that the bar in NGC 4371 was formed at 
a redshift z = been ever since influencing the 

evolution of the galaxy. 

2. We find that the mean stellar age of the more central part of 
the galaxy major disc (our MUSE field does not cover the 
whole galaxy) is above 7 Gyr, with a small contribution of 
younger stars. We interpret this result in the context of the 
evolutionary processes connected to the high-density envi¬ 
ronments in the Virgo cluster. This result thus suggests that 


the removal of gas from the galaxy (ram-pressure stripping) 
and the removal of external cold gas available to fall onto the 
disc and ignite more star formation (starvation/strangulation) 
were already effective at a redshift z = i’ though 
the cluster is still not a dynamically relaxed system today, 
at z = 0. The error bars are derived from the Monte Carlo 
realisations described at the end of Sect. 16.11 

With these data alone, we cannot say much about the fraction 
of barred galaxies at z > 1. However, our results on NGC 4371 
do open up the possibility that - even if bar-driven secular evolu¬ 
tion was not the dominant galaxy evolution process at these red- 
shifts, being overshadowed by galaxy interactions - bar-driven 
secular evolution may have a more extended and significant role 
in the history of disc galaxies, not necessarily restricted to more 
recent cosmic epochs. 

It is evident the necessity of pursuing a similar study as the 
current one, but strengthened with a large, statistically signifi¬ 
cant sample of nearby disc galaxies. Furthermore, in the case of 
NGC 4371 for example, mapping the galaxy to expand the area 
observed with MUSE is potentially very revealing. With such a 
mosaic of MUSE data cubes, the tentative detection of the h 4 
minima along the bar, discussed in Sect. |5.2| could be put on 
firmer grounds, since they seem to occur just at the outskirts of 
our MUSE field. In addition, this work also signals the edge that 
will be gained once the narrow field mode in MUSE is ready, 
providing a spatial resolution of the order of 0.05". At the dis¬ 
tance of NGC 4371, this corresponds to an astonishing spatial 
resolution of only 4 pc, which will clearly unveil the nature of 
the innermost stellar components in the galaxy, and their con¬ 
nected evolution with the central supermassive black hole. 
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